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Cubic H 3 S around 200 GPa: an atomic hydrogen superconductor stabilized by sulfur 

D. A. Papaconstantopoulos^, B. M. Klein^, M. J. Mehl^, and W. E. Pickett^ 

^ Computational Materials Science Center, George Mason University, Fairfax, VA 22030, USA 
^Department of Physics, University of California Davis, Davis, CA 95616, USA 
^Center for Computational Materials Science, Naval Research Laboratory, Washington DC 20375, USA 

(Dated: May 28, 2015) 

The multiple scattering-based theory of Caspar! and Gyorffy for the electron-ion matrix element 
in close packed metals is applied to 7m3m H 3 S, which has been predicted by Duan et al. and 
Bernstein et al to be the stable phase at this stoichiometry around 190 GPa, thus is the leading 
candidate to be the phase observed to superconduct at 190K by Drozdov, Eremets, and Troyan. 

The nearly perfect separation of vibrational modes into those of S and of H character provides a 
simplification that enables identification of contributions of the two atoms separately. The picture 
that arises is basically that of superconducting atomic H stabilized by strong covalent mixing with 
S p and d character. The reported isotope shift is much larger than the theoretical one, suggesting 
there is large anharmonicity in the H vibrations. Given the relative unimportance of sulfur, hydrides 
of lighter atoms at similarly high pressures may also lead to high temperature superconductivity. 


I. BACKGROUND 

The report by Drozdov, Eremets, and Troyan^ (DET) 
of superconductivity up to Tc=190 K in H 2 S compressed 
to the 200 GPa regime has reignited excitement in the 
possibility of achieving room temperature superconduc¬ 
tivity. This report builds on previous success of pressure 
enhancement of T^ in a variety of types of materials: from 
134K to 164K in the cuprate Hg2223,i^ from zero to 20- 
25 K in the simple metals Li, Ca, and Y, and from 
zero to 14K in the insulator siliconJ^ An anticipated ma¬ 
jor factor is the increase in the phonon energy scale with 
compression, since it sets the temperature scale for Tc, 
as pointed out early oni^ and reviewed more recently!^ 
by Ashcroft in predicting possible room temperature su¬ 
perconductivity in metallic hydrogen. 

The newly reported high values of Tc appear to con¬ 
firm theoretical predictions that predated the experi¬ 
ment. Applying particle swarm crystal structure search 
techniques founded on density functional theory, Li et al. 
predicted!^ candidate stable crystal structures of H 2 S up 
to 220 GPa. These predictions were followed by linear re¬ 
sponse calculations of the phonon spectrum uiq , electron- 
phonon matrix elements, and finally the Eliashberg spec¬ 
tral function a^F{ui), from which Tc can be calculated, 
depending only mildly on the presumed value of the re¬ 
tarded Coulomb repulsion /r* = 0.10-0.13. Eor pressures 
of 140-180 GPa, they obtained an electron-phonon cou¬ 
pling strength A=1.0-1.2, an Alien-Dynes characteristic 
frequency!^ Wiog^^lOOOK, and Tc of 40K at 140 GPa and 
peaking at 80 K at 160 GPa. While 80K is well short of 
the reported Tc=190K, the result is convincing that very 
high Tc is predicted in H 2 S at high pressure. 

The sister stoichiometry H 3 S has been explored in 
very similar fashion by Duan et ah^ Predicting struc¬ 
tures to more than 200 GPa, their linear response re¬ 
sults for Im3m H 3 S led to very large calculated values of 
electron-phonon coupling strength A=2.0-2.2, frequency 
scales ujiog ^1300K, and values of Tc up to 200K. In the 
calculations of Li et al. and Duan et al., the high values 


of ujiog are expected from the anticipated increase of force 
constants as volume is decreased. The large values of A, 
a factor of two or more over most other very good super¬ 
conductors (including MgB 2 ), imply that the electronic 
matrix elements are substantially larger than in nearly 
all known conventional superconductors. 

In this report we use Gaspari-Gyorffy (GG) theoryii 
to provide insight into why electron-ion matrix elements 
vary, and evidently increase strongly, with pressure. Such 
understanding is necessary not only to interpret the re¬ 
sults described above, but also to provide essential clues 
how to increase matrix elements, and A, at lower or pos¬ 
sibly ambient pressure. Interestingly, shortly after the 
formulation of this theory, two of the present authors ap¬ 
plied it to predict Tc"^250K in metallic hydrogen at a 
few Mbar pressurej^^ 

GG theoryii builds on the earlier observation of 
Hopfieldi^ that electron scattering off (moving) ions has 
strong local character. Pirst, metallic screening means 
the Thomas-Eermi screening length is short, of the or¬ 
der of 1 A, and very weakly dependent on carrier density 
(more correctly, the Fermi level density of states [DOS] 
N(Ef)). GG employed a multiple scattering Green’s 
function formalism that facilitated three simplifications. 
The first is that the potential is spherical (very good ap¬ 
proximation) and is negligible beyond the atomic sphere; 
the second is that the linear change in potential of a dis¬ 
placed ion can be approximated by a rigid shift of the 
atomic potential. Thirdly, the direction dependence of 
the wavefunction coefficients is averaged out, thereby ne¬ 
glecting any special influence of Fermi surface shape. The 
H 3 S Fermi surface^S is large and multisheeted, minimiz¬ 
ing the likelihood of Fermi surface effects. It is possible 
that these approximations may improve with reduction 
in volume, in any case these approximations should not 
degrade as the system become denser. The bands shown 
by Duan et alJ^ and Bernstein et al^ for the Im3m 
structure of H 3 S that we discuss show much free elec¬ 
tron, spherical character in the lower 75% of the occupied 
bands, though less so around the Fermi energy where S 
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FIG. 1: (Color online) The ImSm crystal structure of H 3 S, 
illustrating the two interlaced ReOa sublattices. Large sphere 
is S, small sphere is H. 


3p - H Is hybridization produces structure in the DOS. 
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FIG. 2: (Color online) Total and orbital-projected densities 
of states of H 3 S at 210 GPa (a=5.6 a.u.). Note the position 
of Ef within the sharp and narrow peak, and that the H Is 
contribution is that from all three H atoms. 


II. THEORY AND RESULTS 


The coupling strength A, and the frequency weighting 
spectrum g{uj) normalized to unity, are given in Eliash- 
berg theory 2 i by 


A = 2 


a^F{uj) 


duj\ 


UJ 


5(w) = 


2a'^F{uj) 
Aw ’ 


( 1 ) 


where a^F is the Eliashberg electron-phonon spectral 
function that governs many superconducting properties. 
The calculations by Li et al. for H 2 S and Duan et al. 
for H 3 S demonstrate that the lower range of phonon fre¬ 
quencies (the acoustic modes) have negligible H charac¬ 
ter, while the optic modes above the gap at 20-25 THz 
have negligible S character, making it an ideal platform 
for applying the GG expression to the atoms separately. 

Thus A = 3Xh + As, where the latter arise from the 
integral over the low frequency S modes, the former from 
the nine higher frequency H branches. In this case the 
GG expression, given originally for an elemental solid, 
can be applied to the S and H spheres separately)^ Each 
atomic {j) coupling constant Xj is given by the integral 
over the appropriate frequency region, leading to 

^ N{Ef) <Ij>_ Pj 
^ ~ M,wf “ M.wf ■ ^ ^ 

The averaged matrix elements </^> obtained from GG 
theory are discussed below. 

The separation of mode character also allows a simple 
estimate of the total frequency moments that enter the 
Alien-Dynes (AD) equation for Tc, through the weight 
function 


5 (w) = [Asgs(w) + 3A_f/g//(w)]/(As + 3Ai/), (3) 

where the partial gj functions, defined analogously to 
that of AD (Eq. [T]) are separately normalized to unity. 


Both H 2 S and H 3 S having been show n^^d^ to have 
strong electron-phonon coupling at high pressure. Bern¬ 
stein et al^ have provided convincing evidence that H 2 S 
is unstable to decomposition into bcc H 3 S and sulfur, and 
that competing stoichiometries are unlikely. This result 
confirms the suggestion of DET, who reported sulfur for¬ 
mation in their samples. Because of this evidence on the 
most likely superconducting phase, we focus on /m3m 
H 3 S with its simple bcc structure based in two interlaced 
ReOs structure sublattices. From a'^F of Duan et al., we 
simplify with a constant a^Fj (constant for each atom 
species j) with frequency ranges (in kelvin) of [430,820] 
for S and [1250,2500] for H. Results are insensitive to 
these limits, depending mostly on the mean frequencies 
and the separation of A into S and H contributions. The 
resulting frequency moments w/og, wi, and W 2 and associ¬ 
ated data, for insight into separate S and H contribution, 
for use in the AD equation^ and to compare with results 
below from GG theory, are presented in Table [U 

With /j,*=0.15, Tc=234 K results; the difference from 
the value quoted by Duan et al. might be due to neglect 
of the strong coupling factor /i=1.13 factor in the AD 
equation, which amounts to a 26 K increase, or partially 
to our constant a^F modeling. Neglecting the contribu¬ 
tion from the S modes, A is decreased from 2.2 to 1.5 but 
ujiog increases from 1500 K to 1770 K, and Tc decreases 
by only 19 K to 215 K. The sulfur contribution to Tc is 
8 %; H 3 S is basically an atomic hydrogen high temper¬ 
ature superconductor. Bernstein et al. also suggested 
that S vibrations are not essential for the high Tc. The 
H isotope effect can also be obtained. A is unaffected by 
masses; Mw^ is a function of the force constants alone, so 
frequencies, specifically Wj^^ decrease as the square root 
of the mass. The resulting critical temperature is re¬ 
duced to ~170 K, slightly more than 234 K/-\/2 because 
the small S contribution remains unchanged. The exper¬ 
imental value of DET— is 90 K; the most likely cause of 
this discrepancy is strong anharmonicity of the H optic 














modes. 

The second moment frequency at 200 GPa of S is 
wf=615 K, while that of H is w|^=1840 K, thus with the 
atomic masses of 32 and 1 a.m.u respectively, the denom¬ 
inator is 32/9 =3.5 larger for S. The consequence is 
that a given contribution to tjh is 3.5x3~10 times more 
effective in increasing A than the same contribution to rjs 
(though in practice there is no clear method of effecting 
such a tradeoff). 


TABLE I: Electron-phonon coupling data for H 3 S, obtained 
from modeling the results of Duan et al. with a constant a^F 
model. Frequencies are in Kelvin. 



S 

H 

H 3 S 

^log 

595 

1770 

1500 

cai 

605 

1800 

1530 

U)2 

615 

1840 

1560 

Mul (eV/A^) 

9.3 

2.6 

- 

d (eV/A^) 

4.7 

1.48 

- 

A 

0.5 

1.7/3 

2.2 


In terms of the phase shifts for the j-th atom for 
orbital channel the square electron-ion matrix element 
averaged over the Fermi surface can be written in the 
simple foriiiii as 

<^j> = ^ NiEp)^ ^ 

where = Nt.{EF)/N^'^ is the ratio of the £-th par¬ 
tial DOS to the single scatterer DOS, for the given 
atomic potential in a homogeneous system. </^> is inde¬ 
pendent of N(Ei^’) since it can equally well be expressed^^ 
in terms of the fractions Ni,{Ep)/N{Ep) which usually 
do not reflect the van Hove singularities of either one. 
The calculated DOS at 210 GPa is shown in Fig. [21 The 
Fermi level falls at a sharp and narrow peak; calcula¬ 
tions at other volumes indicate this is a persistent oc¬ 
currence. <Ij> will tend to be maximized in the cases 
where “neighboring” channels £, £ -|- 1 have a large dif¬ 
ference in phase shifts but ratios Pe that are as large 
as possible. From the GG expression, for each atom 
r]j = N{Ep)<Ij>, and the latter factor involves the 
sin^(5i —factor and products of PDOS ratios pipi+i. 
MjWj can be expressed in terms of the ionic force con¬ 
stants, independent of Mj, (which we return to below) 
so that any isotope effect different from will arise 

from factors beyond A (primarily anharmonicity). 

The calculations have been carried out with two 
all-electron linearized augmented plane wave (LAPW) 
codes, one developed at NRL,— and also ELK.— The 
sphere radii were 1.8 a.u. and 1.0 a.u. for S and H re¬ 
spectively, except for the smallest volume where the S 
radius was reduced because the sum of the radii must be 
no more than a/2. 

The band structur o^^i^° consists of four nearly filled 
bands, leaving some holes at F and electrons around 



FIG. 3: (Color online) Plots demonstrating the very strong 
pressure increase of <I^{j)> (left axis) and rjj (right axis) for 
j=S and H in ImSm H 3 S. For H, the rj plotted is the total 
for 3 H atoms. 


N. In addition, a fifth broad band is roughly half filled. 
The DOS is noteworthy: free electron like over 20 eV of 
the valence band before strong structure arises in a ±7 
eV range centered at EiT-, which lies very close (slightly 
above) the strong and sharp peak in the DOS. This peak 
at Ej? persists for all pressures from P=0 to 300 GPa and 
even above, almost as if E^ were pinned at this peak, 
while other features of the DOS evolve. 

Table El shows the Fermi level values of total and an¬ 
gular momentum components of the electronic densities 
of states across a wide range of volumes. Even though 
/m3m H 3 S may not be stable at lower pressures, we 
provide results for the large range P=0-210 GPa to ob¬ 
serve the effect of interatomic distance on the electronic 
structure and coupling. While the total N(Ei?) shows a 
weak non-monotonic variation, the ^-components have a 
stronger lattice constant (pressure) dependence. As ex¬ 
pected the sulfur 3p-like states are the dominant compo¬ 
nent but a strong participation of 3d character especially 
at high pressures is present, hybridizing with the also 
strong and nearly constant H Is contribution. 

Now we discuss the electronic factor r]j. It should be 
kept in mind that the relative importance for Tc of H 
versus S modes is not simply Aj, but more like XjUJ 2 j 
which is a factor of ten greater for 3H than for S). For 
hydrogen, which dominates the contribution to Tc only 
the s — p channel is relevant at all pressures. The phase 
shift factor sin^((5s — 5p) decreases with pressure. Fig. [3] 
illustrates the factor of two increase in rjH from P=100 
GPa to 300 GPa. This dramatic increase results from 
an even larger increase in the PDOS product , re¬ 

flecting transfer of Is character to 2p character. In the 
spherical harmonic expansion of the atomic wavefunc- 
tions this ‘ 2 p’ character represent the expansion of tails 
of the orbitals on neighboring atoms that gives rise to 
the increased H-S hybridization under pressure. 

As noted above from our analysis of the results of Duan 
et al., the sulfur contribution is less important. As for 










TABLE II: Total and angular-momentum decomposed den¬ 
sities of electronic states at the Fermi level (states/Ry-both 
spins) in ImSm H 3 S. The f-components are projections within 
sphere radii R(S) = 1.8 a.u. and R(H)=1.0 a.u. 


a (a.u.) 

N(Ef) 

S-s S-p S-d 

H-s/atom 

5.6 

6.93 

0.325 1.74 0.751 

0.435 

5.8 

6.43 

0.286 1.53 0.605 

0.405 

6.0 

6.47 

0.274 1.46 0.528 

0.415 

6.2 

6.42 

0.253 1.18 0.528 

0.462 

6.4 

6.79 

0.259 1.36 0.417 

0.455 

6.6 

7.15 

0.255 1.38 0.375 

0.482 

6.8 

7.56 

0.254 1.42 0.339 

0.513 


H, the increase in 77 occurs in spite of a decrease in the 
sin^(5p — 5d) factor, by 30% from P=I00 to P=300 GPa. 
Over this pressure range, the PDOS ratio product 
increases by 75%, giving a net rj increase by more than 
50%. In addition, the S s—p and even d—f channels begin 
to contribute, reaching 20% of the total of <I^> at 300 
GPa. Thus the increase is a composite effect of increase of 
both d and s character of S, which is reflected also in the 
growing phase shifts of these channels. These transfers 
of atomic character under pressure are consistent with 
general expectations of the evolution of atomic character 
under reduction of volume. Interpolating to 200 GPa to 
obtain 775 = 5 .84 eV/A^, 7757 =!.71 eV/A^ (which must be 
multiplied by three), and using the frequency moments 
from Table I, we obtain A=2.6, thus Tc 250K depend¬ 
ing somewhat on the chosen value of 77 *= 0 . 10 - 0 .I 3 . This 
value is very consistent with the value of Tc (above) ob¬ 
tained from our modeling of the Duan data. Though the 
numbers might vary for other low-Z hydrides (H„P, H„B, 
...) the lack of any special role of S in these results sug¬ 
gests there should also be a strong increase in Tc with 
pressure in other low-Z element hydrides. 

A question of great interest is whether Tc increases fur¬ 
ther at higher pressure. The data presented in Table Hill 
provides the pressure dependence of the important quan¬ 
tities entering 775 and rjH- This data demonstrates that 
the strongest contributions arise for sulfur from the pd 
channel and for hydrogen from the sp channel. Note 
that for H one should multiply by three to account for 
the three H atoms in the unit cell. 

Since Xh dominates the sulfur contribution, we can 
focus on the H contribution alone. The total pressure 
derivative contains several contributions 

dlogTc _ dlogojiog dlog/i dlogA(A) 
dP ~ dP ^ dP ^ dP 

_ dloguJiog , fdlogfi ^ d\ogE{\)^d\ 

~ + dx ’dP' 

where /i is the strong coupling correction and E{X) is 
the exponential term in the Alien-Dynes equation. The 
pressure variation of the Mujfp denominator is challeng¬ 
ing to approximate without full calculations of the spec¬ 
trum and a^E. We have modeled the variation of the 
H spectrum by assuming the three 3-fold P point optic 


modes are representative. The lower two of these modes 
are IR-active involving H-S bond stretch and bond bend¬ 
ing modes, the hardest frequency is a silent mode with 
quadrupolar H motion with respect to S. 

From calculations of these P frequencies in the 240- 
270 GPa range using the ELK code we calculate a posi¬ 
tive but modest pressure increase dlogun/dP~l.'dxlQ~^ 
GPa“^. From Fig. we obtain dlogrjH/dP=‘ii.?>x\Q~^ 
GPa“^, thus A decreases with pressure approximately as 
dlogA/ 7 /dP~- 0 . 3 xI 0 “^/GPa. However, to our precision 
this is indistinguishable from zero, so the pressure deriva¬ 
tive in Eq. [5] reduces to the first term, the frequency 
derivative. The resulting prediction is a small increase 
dTc/dP=QA K/GPa. This result disagrees in sign with 
Duan et al., who quoted a smaller (in magnitude) nega¬ 
tive value of dTc/dP=-0.12 K/GPa from direct calcula¬ 
tion, however both numbers are small compared to the 
large value of Tc itself, so there is no significant disagree¬ 
ment. 


III. CONCLUSIONS 

The report by Drozdov, Eremets, and Troyan of Tc 
up to 190 K in H„S samples has breathed new life into 
the 50 year old expectation of high Tc in atomic H sys¬ 
tems. Both Li et al. and Duan et al. had found that 
Eliashberg theory and linear response results for electron- 
phonon coupling account for Tc in the 80-200 K range 
for H 2 S and H 3 S at high pressure, and the analysis of 
Bernstein et al make ImAm H 3 S the primary candidate 
to be this record-high temperature superconductor. In 
this paper we have established that the coupling of H vi¬ 
brations increases strongly for pressures up to and even 
beyond 210 GPa, and that 90-|-% of the coupling arises 
from H vibrations in this hydride that is stabilized by 
hybridization^^ with S. This picture is analogous to the 
finding of the essential contribution of H in the super¬ 
conductor PdH at ambient pressure and the broader 
picture of Ashcrofti^ of superconducting elemental H at 
high pressure. The theoretical isotope shift of Tc based 
on the harmonic approximation is not in agreement with 
the experimental result, suggesting substantial H anhar- 
monicity will be necessary to understand before the pic¬ 
ture is complete. Our picture, which relies on coupling 
across the large Fermi surface, is at odds with the hole 
superconductivity picture of Hirsch and Marsiglioi^ 
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TABLE III: Pressure variation of the Hopfield parameter 77 (eV/A^) and electron-ion squared matrix element <7^> (eV^/A^), 
both decomposed into the three channels sp, pd, and df for S, and the H sp channel (given for a single atom). 


a (a.u.) 

P (GPa) 

S Tjsp 

S Rp 

S f]pd 

SlU 

S Pdf 

Sidf 

H T}sp 

Hlip 

5.6 

210 

0.68 

2.66 

5.12 

20.1 

0.31 

1.21 

1.80 

7.06 

5.8 

142 

0.39 

1.65 

4.45 

18.9 

0.21 

0.89 

1.43 

6.06 

6.0 

92 

0.24 

1.00 

4.14 

17.4 

0.15 

0.63 

1.24 

5.21 

6.2 

57 

0.19 

0.78 

4.13 

17.5 

0.13 

0.57 

1.02 

4.34 

6.4 

31 

0.08 

0.33 

3.52 

14.1 

0.08 

0.31 

0.96 

3.87 

6.6 

13 

0.05 

0.18 

3.37 

12.7 

0.06 

0.21 

0.87 

3.31 

6.8 

0 

0.02 

0.09 

3.12 

11.2 

0.04 

0.14 

0.79 

2.83 
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